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Abstract 

We consider a nonhomogeneous Burgers equation with time variable coefficients of the 
form Ut + {ij.{t) / fj,{t))U + UUx = {l/2fj,{t))Uxx — uj'^{t)x, and obtain an explicit solution of 
the general initial value problem in terms of solution to a corresponding linear ODE. Special 
exact solutions such as generalized shock and multi-shock solitary waves, triangular wave, 
Af-wave and rational type solutions are found and discussed. As exactly solvable models, 
we study forced Burgers equations with constant damping and an exponentially decaying 
diffusion coefficient. Different type of exact solutions are obtained for the critical, over and 
under damping cases, and their behavior is illustrated explicitly. In particular, the existence 
of inelastic type of collisions is observed by constructing multi-shock solitary wave solutions, 
and for the rational type solutions the motion of the pole singularities is described. 

1 Introduction 

The nonlinear diffusion equation, known as Burgers equation (BE) after the extensive work of J. 
M. Burgers, [3 [5] is an important model which appears in various fields of physical science. In 
hydrodynamics, it is a standard model of turbulence used to study propagation of nonlinear waves 
and shock formation, [3 . It is used also to describe processes in gas dynamics [U [S] , nonlinear 
acoustics [6], heat conduction, and plasma physics. In cosmology, the Burgers equation is a good 
approximation to understand the formation and distribution of matter at large scales, [7]. 

The standard Burgers equation is of the form Vt + VVx = vVxx, where V mostly represents 
the velocity field, t is a time variable, a; S M is the space variable and z/ is a constant viscosity 
or diffusion coefficient. This equation is probably the simplest nonlinear model admitting direct 
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linearization, and thus being C-integrable in contrast to S-integrable systems which require spectral 
transform technics. Indeed, exact explicit solutions of the Burgers equation can be obtained by 
the Cole-Hopf transformation, which transforms the nonlinear Burgers equation to a linear heat 
equation, [4j[5]. Beside that, simple Ansatz method was also applied to find special solutions, like 
traveling waves and similarity solutions. Lately, other methods like, Hirota's direct method and 
Backlund transformation hyperbolic function method 8 , and homogeneous balance method 
[5] were used to construct new exact solutions of the BE. 

As known, the main features of the Burgers equation are due to the simultaneous existence of 
a nonlinear term and a linear diffusion term. If the diffusion is dominant over nonlinearity, the 
solution of the BE approaches the solution of the diffusion equation. On the other hand, if the 
nonlinear term dominates over the diffusion, one may expect formation of shock discontinuities. 
An interesting property of the BE appears when a balance occurs between the nonlinear effect 
and the effects of dissipative nature. In that case, the system exhibits shock profile solitary wave 
solutions. Moreover, it is known that Burgers equation has also multi-shock solitary wave solutions 
[TOl[Tl]. In that case, shocks of different amplitude and speed can fuse (merge) to a single shock, 
so that completely non-elastic interactions may occur. Another important property of the Burgers 
equation is related with the rational type solutions. It is well known that, the zeros of the heat 
equation solution lead to pole singularities for the Burgers solution. Choodnovsky brothers |12j 
and F. Calogero .13,, showed that the motion of the poles corresponds formally to the motion of 
one-dimensional particles interacting via simple two-body potentials, such that the corresponding 
many body problems are integrable. For recent work on the pole dynamics of the standard Burgers 
equation one can see |14j . 

The standard BE, as mentioned above, is a well known exactly solvable model. However, the 
nonhomogeneous and variable parametric Burgers equation, in general, is not integrable and very 
few exactly solvable models are known. For example, in case Ut + UUx = vUxx + F{x,t), if the 
forcing term depends only on time, i.e. F(x,t) = G{t), this equation can be transformed to a 
standard Burgers equation, see |T5]. The IVP with an elastic forcing term F{x,t) — ~k^x -f 
f{t) is discussed and analytic solutions are obtained in [TB]. Later, the problem with F{x,t) = 
G(t)x, where G{t) is arbitrary, was completely solved in terms of solution to the standard Burgers 
equation, see [171 . In [18| , an invertible transformation between the nonhomogeneous BE and the 
stationary Schrodinger equation was constructed so that each solution of the stationary Schrodinger 
equation generated a fully time-dependent solution of the nonhomogeneous BE. Recently, exact 
solutions were obtained using the Cole-Hopf transformation and the Green's function approach, 
see [19] . Transformation properties of a variable-coefficient Burgers equation were discussed in 
[2U] . In PJ], a forced Burgers model with space- and time-dependent coefficients of the form 
Ut + a(x, t)UUx = b{x, t)Uxx + F{x, t), was investigated using a generalized Cole-Hopf transform 
and symbolic computation. For the significance of the generalized forced Burgers models and 
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possible applications in various fields one can see again the discussion in '2T and references given 
there. 

In this work, we consider a nonhomogeneous Burgers equation (NHBE) with time variable 
coefficients of the form 

Ut + + UU^ ^ -^U^^ - uj'^{t)x , -cx)< a; < CX3 (1) 

where T{t) — fi{t)/ii{t) is the damping term, D{t) ~ l/2fi{t) is the diffusion coefficient, and 
the forcing term which is linear in the space variable x. In Sec. 2, we show 
that solutions of the NHBE ([T} can be obtained in terms of solutions to the standard BE or heat 
equation and a related linear ODE. As a result, an explicit solution for the IVP of the NHBE 
([T|) is found in terms of solution to the corresponding second order linear ODE with variable 
frequency and damping. Then, some particular exact solutions such as shock and multi-shock 
solitary type waves, triangular wave, N-wave and rational type solutions are obtained. In Sec. 3, 
for comparative reasons, first we recall some solutions of the nonhomogeneous Burgers equation 
with constant coefficients. Then, exactly solvable NHBE models ^ with positive constant damping 
and exponentially decaying diffusion coefficient are considered. Different type of exact solutions 
mentioned in Sec. 2 are obtained for the critical, under and over damping cases. We observe 
generalized traveling wave solutions which speed, steepness, and shock amplitude are functions of 
time. Special properties like interaction of shocks in multi-shock solitary solutions, and motion of 
pole singularities of rational type solutions are described explicitly. Sec. 4 includes brief summary 
and future plans. 



2 Nonhomogeneous Burgers equation with time dependent 
coefficients 

In the following proposition we obtain relation between solutions of the nonhomogeneous variable 
coefficients Burgers equation and the standard Burgers equation. Then, Cole-Hopf transform 
allows us to find an explicit solution of the IVP for the variable coefficient NHBE ^ in terms of 
solution to a corresponding second order linear ODE. 

Proposition 2.1 // r{t) is solution of the IVP for the linear ODE 

r+|^r + c^'(t)r = 0, r(<o) - ro ^ 0, r(io) = 0, (2) 

then the IVP for the NHBE with variable coefficients 



U{x,t)\t^to =U{x,tQ) , -oo < a; < oo 
has solution in the following forms: 



(3) 



a) U{x,t)='-^x + ^^V{^{x,t)Mt)). (4) 

r{t) ti{t)r{t) 
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where 



r{t) ' '"'J f^iOrHO 

and the function V{ri,T) satisfies the IVP for the standard BE 

F(r,,0) =/x(to);7(7?,to). 



T{to) = 0, 



(5) 



(6) 



b) 



(7) 



rit) ^i{t)r{t) ip{T]{x,t),T{t)) ' 
where 77, r are as defined in part (a), and (p{r},T) satisfies the IVP for the heat equation 

^(77,0) =exp [-pfiito)U{^,to)d^]. 

Proof: a) Using the Ansatz U{x,t) ^ [fi{t)]-^[p{t)x + s{t)V{s{t)x,T{t))], it is easy to show 
that, if the auxihary functions satisfy the nonhnear system of ordinary differential equations 



(8) 



fi{t) 



T- -777=0, T(io)=0, 



(9) 



H{t) 

then the IVP Q for the NHBE transforms to the IVP ^ for the standard BE. Also, noticing that 
Eq.® is a nonlinear Riccati equation, the system is easily solved, and we obtain the functions 



p{t) - fi{t) 



■it) 



^{t)^r^to) 



, T(to) = 0; 



sit) 



r{to) 



(10) 



which substituted back in the Ansatz give the result Thus, solution of the NHBE ^ is 

explicitly obtained in terms of solution V{r], r) to the BE ([5]) and solution r{t) of the IVP for the 
linear ODE ([2]). Part b) of the proposition, follows directly from the Cole-Hopf transformation 
V = —(pjj/ip, which reduces the IVP (jH)) for the BE to the IVP © for the usual heat equation. □ 
As well known, the IVP ([5]) for the heat equation has solution 



(p{r),T) 



/2ttt J- 



exp 



2t 



^(e,o)de, 



and Colc-Hopf transformation V ~ —ifri/f, leads to solution of the IVP ([6]) for the BE 
V{v,r) = 





exp 


- 


[^ + fv{C,o)de) 


' d^ 


rOG 

Loo 


- {^ + l'vie,o)de) 


d^ 



Therefore, using the above proposition, one can find formal solution of the IVP ([3]) for the NHBE 
in terms of solution r(t) of the linear ODE that is 



Uix,t)^^-^x- 



r{to) 1 






1 exp 


- 


{j^^^^ff^ + f Kto)U{C,to)dC^ 




d^ 


p{t)r{t)\ 


rOO 

J-oo exp 


- 


{j^^ + ft^ito)U{i',to)d^'^ 




di 



,(11) 
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where T{t) is as defined in ([5]), and the time interval on which the sohition exists depends on the 
properties of the auxiliary functions. Since it is difficult to analyze solution for an arbitrary 
initial condition, in what follows, we consider particular problems for which the NHBE ^ subject 
to some localized initial profiles has exact solutions and one can observe explicitly their behavior. 
As known, the standard BE (jH]) has different type of solutions, such as shock solitary waves, 
similarity, N-wave and rational function solutions. This suggests us to look for corresponding type 
of solutions for the NHBE with variable coefficients, as follows. 

a. Shock solitary wave solutions. The standard BE ([6]) has shock solitary wave solution 

V{t],t) =c- Ata,nh[A{r]-cT + co)], (12) 

where A, c, cq are arbitrary real constants. Then, the NHBE ([3]) with initial condition 

U{x, to) — ^ (c — Atanh [Ax + cq]) 

has generalized shock solitary wave solution 

" 'rito] 



U(x,t) = 4^2^+ -7T^ ( c-^tanh 
^ ^ r(t) f,{t)r{t) * 



X — CT{t) + Co 



(13) 



T{t) 

where r{t) is a solution of the IVP ^ and T{t) is given in ([SJ. In particular, when in ([T^ one has 
c = Co — 0, the NHBE ^ has shock type static solution of the form 

'Ar{to) 



^ r{t) Ar(to) 

U(x, t) = -^x --^^ tanh 

^ ^ r{t) fi{t)r{t) 



(14) 



r{t) 

corresponding to the initial condition U{x,to) — —{A/ fi{to)) tanhlAx] , — oo < a; < oo. 

Solution (IT^ of the standard BE © is a localized wave of constant amplitude moving with 
constant speed. However, the solution (IT^ of the NHBE with variable coefficients is a generalized 
traveling wave of the form U{x,t) = u{x,t) + A{t) tanh[B{t){x — C{t))], where the term u{x,t) 
contributes to the wave amplitude, A{t) is the shock amplitude, B(t) is related with the steepness 
of the shock profile, x — C{t) describes the motion of the "center" of the profile, and v — C{t) 
is its velocity. Accordingly, for the wave solution (IT51) . one can see that the shock amplitude is 
proportional to l//i(i)r(t), and steepness of the profile is proportional to \/r{t). Also, the position 
of the "center" of the wave profile is described by x{t) — [r{t) / r{to)\{cT{t) — cq), where its velocity 
can be easily found using that v{t) — x{t). 

Multi-shock solitary wave solutions. Since the standard BE has multi-shock solitary wave 
solutions, [ini E], it is natural to ask for multi-shock solitary type solutions for the NHBE with 
variable coefficients. Here, we outline the procedure, and give formal results. Clearly, the heat 
equation has simple solutions of the form 

^2 

fi{v,T) = exp[pi{T],T)] , p^{T],T)^-atr]+Y'^+Pi, aj,p- G K, 
and their linear superposition 

(p{ri,T) = exp[pi(77,T)] +exp[p2{v,T)] + ... + exp[pk{T],T)], i = 1,2, ...,fc. 
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is also a solution. By the Cole-Hopl transform V = 
corresponding solutions of the form 

aiexp[pi(?7, t)] + 02 exp[p2(?7, t)] + 



1/(77, r) 



(P',j/ip, it follows that the BE ([6]), has 
-afcexp[pfc(?7, r)] 



exp[pi(7],r)] +exp[p2{r],T)] + ... + exp[pfe(?7, r)] 
Therefore, using Proposition 2.1, one obtains that the NHBE ([3]), has generalized multi-shock 
solitary wave solutions given by 



where 



ai exp[pi(a;, t)] +02 exp[p2(a;, i)] 



afcexp[pfe(x,t)] 



exp[pi(a;,i)] + exp[p2(a;, t)] + ... + cxp[pfc(a:, t)] 



i — 1, 2, fc. 



r(t) is solution of ([2]), and r(i) is given in ([S|). For k — 2 and proper choice of constants, one 
obtains one-shock solitary wave. When fc > 2 one expects formation of multi-shock solitary type 
solutions. Indeed, this is the case, and illustrative examples are given in Sec. 3. 

b. Triangular wave solution. The BE © has triangular wave (similarity) solution, 

1 / (e^ - 1) exp[-7yV2r] \ 



y(77,r) 



(15) 



/2^ \^ 1 + i(e^ - l)erfc[T]/V2^] J ' 
corresponding to initial condition V{ri, 0) — A6{ri), see [TU], where A is a constant, 6{ri) is the Dirac- 
delta distribution, and erfc[a] ~ /^exp[— . Then, the NHBE ([3]) has generalized 

triangular wave solution of the form 

{eA — 1) exp 



U{x,t) = 



rito) 



-(^^)V2r(t) 



r{t) ^i{t)r{t) ^2TTT{t) \^ 1 + 1 (eA _ i)^^^^ VM*)" 
c. iV-wave solution. The heat equation ([5]) has solution 

(p{r],T) = 1 + \-exp[~7f/2T], T>0, 



-{t) > 0. (16) 



which behaves like delta distribution as r — > 0, (a-positive constant). The corresponding N-wa^ve 
solution of BE see [TO], is 



V{v,t)=(1 



r > 0, 



y^a/r exp[— ?7^/2t] 
-tJ I + yfajr exp[— 772/2T] ' 

Therefore, generalized A^-wave solution of the NHBE ^ is of the form 



T{t) > 0. (17) 



r{t) ii{t)T{t)r-^{t) 1 1 + y^77(I) exp[~(^a;)72T(t)] 

Since the behavior of this solution at T(to) = is rather complicated, as an initial profile one can 
consider a profile at any time t > Iq. 

d. Rational function solutions. Formal solution of the IVP for the heat equation ([5]) can 

be found also by applying the evolution operator to the initial condition, that is 

00 

^(r;,r) =exp[(r/2)a2]^(^,0) =^-(^)'=9f ^(r;,0). 



A:=0 
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If the initial condition is (pmiVj 0) = V™i then the solution of the heat problem ([8]) is 

^™(??,r) =exp[(T/2)92],,™ ^iJ„(^,r/2), m = 0,l,2,..., 
where Hm{'ri,T/2) are Kampe de Feriet polynomials, defined by 

with [m/2] = to/2 for even m, and [to/2] = (to — l)/2 for odd to, see j^S]. Using also the relation 
d,^H,n{ri,T/2) ~ mH„i-i{r],T/2), it follows that, the BE ^ has rational solutions of the form 

V„i{r],T) = -—[llKpmiV,^)] = ^ ^1^^ , TO =1,2,3,..., (18) 

or more generally 

/ N ^ (: f \\ T,t=i ma^H^.i{r],T/2) 
Vk{v,r)^-Tr[^^2^a^V^{'n,T)\^ J^—^ fc = 1,2,3,..., (19) 

where Um are arbitrary real constants. Therefore, we obtain that the variable coefficient NHBE 
([3]) with initial conditions Um{x^t{)) = — (TO//i(to)x), to — 1,2,3, has rational solutions 

and with more general initial conditions 



it has rational solutions of the form 

The solutions (P(7| and (PT|) can be written also in terms of standard Hcrmite polynomials Hm{y) 
using that 

(,/t/2)"^ ( in \ 
H,n{v,r/2) = / ^ H,n , i/™(?7,0) =77", to = 0,1,2,..., (22) 



^2^t/2^ 

where Hm{y) are defined by exp[2?/^ — = X]m=o(^"V™')-^m(y)- Thus, the points where Kampe 
de Feriet polynomials vanish can be found in terms of the well known zeros of the Hermite poly- 
nomials. For this, we denote by , / = 1,2, ...,ni the zeros of the Hermite polynomial Hm{y), 
so that for each fixed to, one has HmiVm) = for all I = 1,2, to. From relation it follows 
that 

= -7 97/(' 



iy™(77,T/2) = ^ ii = -i2yli>^T/2, / = 1,2,...,to. (23) 
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Thus, Um{x,t) given by ((20|) has singularities at points where Hm ((r{to)/r{t))x,T(t)/2) ~ 0, and 
according to , the motion of these pole singularities is described by 

.«(0 = -^2y«^y^, Z = l,2,...,m. (24) 
r(ta) V 2 

We note that, for a real- valued solution r{t) and T(i) > 0, t E I, the solution Um{x,t) does not 
have moving singularities on the real line. It may have real singularity only at a; = 0. On the other 
hand, for some special choice of the coefficients a,„, the solutions Uk{x,t) of the form (PT|) may 
have singularities moving on the real line. Illustrative examples are given in next section. 

3 Exactly solvable Burgers models 

3.1 Forced Burgers equation with constant coefficients 

Burgers equation with constant coefficients and a forcing term linear in the space variable x 

Ut + UU, = ^U^^- Lolx, (25) 

is a known integrable model and one can see for example |17[ 124) . For this model, according to 
Proposition 2.1 one has /i(t) = 1,7 = 0, and w^(t) = cjg — real constant, so that the corresponding 
IVP for the second order linear ODE is 

f{t) + ujlr{t) = 0, r(0) = ro 7^ 0, f(0) = 0. (26) 

When, cjq = 0, one has r{t) — j-q, t] — x, t = t, and formula (H)) gives U{x,t) = V{x,t), which is 
a solution of the standard Burgers equation, as expected. Using the approach in previous section, 
we will recall some particular solutions for the cases when ojq > and luq < 0. 

3.1.1 Case uj^ > 0. 

In that case the IVP ((25)) has oscillating solution r{t) = tq cos(a;oi) and the auxiliary function is 
T{t) = ta.n{bJot)/ujo. From Proposition 2.1, it follows that the forced BE ([25]) has solutions 

U{x,t) = —LUotaii{ujot)x + sec{ujot)V{sec{ujot)x,taii{ujot)/ujQ), (27) 

where V{ri, r) is a solution of the standard Burgers equation. In what follows, using the discussion 
in previous section, we will write explicitly some special solutions of BE (j25p . and note that, in the 
limit case — > 0, these solutions U{x,t) approach the solutions V{x,t) of the standard BE. 

a. Forced Burgers equation ([25]) with initial condition U{x, 0) = —Atanh{Ax) , — oo < a; < oo, 
has shock type static wave solution 

U{x, t) — —ojQ ia,n{ujQt)x — A sec(a;oi) tanh {A sec{ijOot)x) , (28) 
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and with initial condition U{x^ 0) — c ~ Atauh[Ax + cq] , oo < a; < oo, it has shock sohtary type 
solution of the form 



(29) 



U{x,t) = —ujQ ta,n[ujQt)x + sec(a;ot) c — A tanh 



A sec{uJot)x tan(cjo<) ) + co 



b. Forced BE ([25)) has triangular wave solution 

/ 



U — — wq tan(a;oi)a; + sec(ti;o^) 



27rtan(ajot) 



{e^ — 1) exp 



2 tan(a;oi) / cjq 



V 



1) er/c 



for t G (0, 7r/2a;o). In general, triangular wave exists on time interval where T{t) — tan(ajo^)/'^o > 0- 

c. A^-wave solution can also be considered for BE ([^5)1 using formula ([T7)) on time interval 
where r(t) > 0. 

d. Forced BE (1251) subject to initial conditions Um.{x, 0) — —m/x , m — 1,2,3, has rational 
function solutions of the form 

mHm^i {sec{ujQt)x,tan{ujQt)/2LUQ) 



Hm {sec{uJot)x, tan(a;oi)/2aJo) 



Um{x,t) = ~ujQtan{LUQt)x — 8ec{ujot) 
with moving pole singularities Xm{t) = —iy\n cos(ujQt) ^ 2 tan(a;oi) / ^^o , ^ = Ij 2, m. 
3.1.2 Case uJ^ < 0. 

When = -i^"^, uj > 0, the BE becomes 
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Ut + UUx^ -Uxx+<^ X, 



(30) 



and the corresponding ODE (j26p has solution r{t) — tq cosh((I't) , where the auxiliary function is 
T{t) ~ ta.nh{ojt) / UJ . Therefore, solutions of BE ([50]) are of the form 

1 , , / X tanh(cDi) ' 



U {x, t) — Lj tanh((I)i)x + 



-V 



(31) 



cosh(wi) \cosli(wi)' uj 

where V{'q,T) is a solution of the standard BE. As in the previous case, we see that in the limit 
a) — ?> 0, the following particular exact solutions U{x, t) approach the corresponding solution V{x, t) 
of the standard BE. 

a. The forced BE ([50)1 with initial condition U{x,{)) ~ c— Atanh[yla:], —co < x < oo has 
one-shock solitary type solution 



U{x, t) = u) ts.n\\(u)t)x + 



1 



c — A tanh 



A 



cosh(a;i) 



sinh((I)<) 
[x-c ) 



cosh(a)i) 

which amplitude depends on time, the center of the wave profile moves according to x[t) 
csmh.{Cjt)/uj, and its velocity is v{t) = ccosh(wi). 

In general, BE pOI) has multi-shock solitary type solutions of the form 



(32) 



U{x, t) = Cj tanh(wt)a; + 



1 



cosh((I;t) 



ii exp[pi(a;,t)] + a2 exp[p2(a;, t)] 



afcexp[pfe(x,i)] 



e-xp[pi{x,t)] -I- exp[p2(a;,i)] + ■•• + exp[pfe(x, t)] 



, (33) 
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where 



Pi{x,t) 



cosh(a)i) 2a) 
b. Triangular wave solution of BE ([50)1 is 



a; + ■^tanh(a)i) i = l,2, ... 



U{x,t) — w tanh((u)i)a; ■ 



cosli(a)i) y 27r tanh(wi) 



V 



(e"^ — 1) exp 


f X "(2 // 2 tanh(wt) \ 
Vcosh(it) / V li 




1 + ^{eA - l)er/c 


/ K W / 2 tanh(wt) 



which corresponds to Dirac-delta initial distribution, 
c. N-wave solution of BE ([50)1 is of the form 



U {x, t) — Co tanh(a)<)x + 



sinh((i)t) cosh((ii<) 



tanh(a;f 



exp 



( X \2 /tj/ tanh(i>£)^N 



1 



■ exp 



'(cosh(it))^/2( 



taiih(;I)i) 



(34) 



tanh(cjf) 

d. Rational function solutions of BE ([5(1)) subject to initial conditions Um{x, 0) = —m/x , m 
1, 2, 3, are of the form 



Um{x, t) = Co taTi]i{Cjt)x + 



1 



mHm-i {xj cosh(ci)i), tanh((I)t)/2;D) 



(35) 



cosh(a)i) \ i7„i (a;/ cosh(a)t), tanh(wf)/2cj) 

For each fixed m, the rational solution (j35p has pole singularities whose motion in the complex 
plane is described by 



a^^^W = -*y£^cosh((I.i)A/2(^^^^i$!^), t>0, ? = l,2,...,m, 



and as (2; they approach the well known poles Xm{t) ~ —iym\j2t of the standard BE. 

Finally, we note that for the above particular solutions, one has limt_).oo V (x, t) ~ Cox. Similar 
result was obtained in [22] and [23] , where the long-time asymptotics for solutions of the BE (|30|) 
were discussed according to the properties of the initial profile. 
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3.2 Forced Burgers equation with constant damping and exponentially 
decaying diffusion coefficient 

In this part, we consider exactly solvable forced Burgers equations of the form 

Ut+-fU + UU^ = ]^e-'''^U^^-oolx , -CX3 < a; < oo, i > 0, (36) 

with constant damping T{t) = 7 > 0, exponentially decaying diffusion coefhcient D{t) — e"''^ /2, 
0J^{t) — cOq > and ii{t) = e''*'. The corresponding IVP for the linear ODE is then 

f + 7f + cj^r = 0, r(0) = ro 7^ 0, f(0) = 0, (37) 

and it has three different type of solutions depending on ujq — (7^/4). In what follows, for each 
case we discuss separately the related variable coefficient Burgers equations ((36]) . 
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3.2.1 Critical damping case. 

If ujI - (7V4) = 0, the IVP dSZl) for the linear ODE, has solution 



ri(0=ri(0)e-*(l + |i), 



(38) 



and the auxiliary function is T\(t) = +7i/2). Therefore, the BE ((36)) has solutions of the form 



TTft 



V 



(39) 



where V{'ri, r) satisfies the standard BE. 

Clearly, when ujq — > 0, one has also 7 — > 0, and in that case we can see that the following 
particular solutions U{x,t) approach solutions V{x,t) of the standard BE. 

a. BE with initial condition U{x,0) = c — Atanh [Ax] has shock solitary type solution 

t \ /e-^*/2- 



c ~ A tanh 



^(^(-c.e--/^) 



(40) 



which shock amplitude decays with time eventually going to zero, and its "center" moves with 
velocity v{t) = c(l - 7t/2)e~''*/^, see Figffl Muhi-shock solitary type solutions of the BE 




Figure 1: Critical damping case. One-shock solitary type solution U{x,t) given by (|40)) . j — 2, 
A = 4, c = 5,co = 0. 



can be found from the general solution 



i + it 



where 



ai e^p[pi{x,t)] 



akexp[pk{x,t)] 



exp[pi{x,t)] + ... + exj>[pk{x,t)] 



(41) 



i + ii 



l,2,...,fc, 



and ai,p^ are real constants. In Figl2]we plot two-shock solitary type wave solution U{x,t), with 
special choices in (|H|), fc = 3, 7 = 2, ai = 1, p° = 0, 02 = 7, = 5, 03 = 15, p§ ~ —A. We observe 
fusion of two-shock solitary wave, which shock contribution eventually goes to zero. 
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0.0 



Figure 2: Critical damping case. Fusion of two-shock solitary type solution U{x,t). 

b. Using again the results in Sec. 2, one obtains that the BE has generalized triangular 
wave solution 



t \ e~i' /(T+ jt/2) ( je^ - 1) exp [-e^*xV2t(l + -/t/2)] 



2 V(l + i<)/ (l + ii)V 27Tt i(pA-l)er/c[e')'*/2a;/y2i(rT7i72) 

corresponding to Dirac-delta initial profile. 

c. BE ((36)l for the critical damping case has A^-wave solution of the form ((T7| with r{t) = ri(t) 
and T{t) — Ti{t). Explicit form of this generalized A^-wave solution is 



2 \{l+-jt)J til + ii) \^ 1 + y^O^exp [-e^*xV2t(l + jt/2)] 

d. Rational solutions of the BE (l36|) with initial conditions Um{x, 0) = —m{l/x), m — 1, 2, ... 
are of the form 

According to p4l) . for each m, the motion of the pole singularities is described by 



xW(i)=-*y2 2;We-^*/V^(l+7V2), ^ = 1,2,.... (43) 

Since 7 > and t > 0, clearly Um{x,t) has no moving singularities on the real line. On the other 
hand, forced BE p6p with more general initial conditions 

k k 
m— 1 m— 

has the fohowing rational solutions 

UJ. t) - -{I f ( ' ] X - ( ^^-i"^°-^--^(t^^' whr), , 

^ ^ ^ \ Z^m=0 (^Th^^' 2(l+it) 
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At the end of this section, we note that for the above particular sohitions one has Imit^oo U (a;, t) = 
— {^/2)x, so that in the fong-time hmit the system becomes stable with velocity proportional to 
the displacement. 



3.2.2 Under damping case. 

If - (7^/4) > 0, then the ODE has oscihating solution 

r2it) = r2(0)^^e-^ cos[nt - a], 



where fl — ^ljq ~ (7^/4) > 0, a = tan ^i-^), and the auxiliary function is 



UJn 



(45) 



(46) 



Then, BE (|36p has solutions of the form 



+ - 



-7t/2 / Q 



=7t/2 



(tan[r!t~a] + ^) 



(47) 



loq cos[f2i — a] \loo cos [fit ~ a]' loq 
where V{ri, r) satisfies the standard BE. 

When 7 — > 0, one has — > cjq > and a — > 0. In that case, it is not difficult to see that the 
bellow given solutions of the forced BE (j36l) with variable coefficients approach the corresponding 
solutions of the forced BE ([25]) with constant coefficients, 
a. BE ([55)) has shock type static solution 



C/(x, ^) = - (I + 1^1 tan[m -a]^x 



An 



cjQ cos[f2t — a] 
and it has shock solitary type solution in the form 



tanh 



U{x, t) ^ - i^^+n tan[nt -a])x 



An 



cos[r2t — a\ 



(48) 



n 



ujf) cos[r2< — a] 



A tanh 



n 



alt/2. 



n 



Wo COS [fit — a] 



tan [fit — a] 



2nJ 



(49) 



Note that the solitary wave is broken by shock discontinuities which appear periodically at finite 
times. Multi-shock solitary waves can be obtained using the solution 

aiexp[pi(x,t)] +02 exp[p2(2;,i)] + ... + akexp[pk{x,t)] 



Uix,t)^t4), + e~^tr_2{tol 



r2{t) 



r2{t) 



where 



Pi{x,t) 



exp[pi(a;,t)] + exp[p2(a;, t)] + ... + cxp[pfc(a:^, 0] 

i — 1, 2, fc, 



r2{t) is given by pS)) . and T2{t) is given by (pS)) . 

b. BE ([55]) has generalized triangular wave solution of the form ([T5| with r{t) — r2{t) and 
r(t) — T2{t) on a time interval where T2{t) > 0. 

c. Generalized A^-wave solution of BE (p6)) for the under damping case can be found using 
formula (|17p on a time interval where T2{t) > 0. 
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d. BE (1361) with initial conditions 

k 



Uk{x, 0) = -(^ ma^x"-!)/ J2 amx"", k ^ 1,2, .. 

m—1 m—O 

has rational function solutions 



3.2.3 Over damping case. 

When uj^ - (7V4) < 0, the IVP ^ has solution 



^3(i)=^3(0)^e-^sinh[r!'i + /3], (51) 



where O' = ^[wg - and ^ = coth"'(2^). Then, 



^'^^^ = "5 (coth[17'i + /3] - ^) , (52) 
and thus BE (^S)) has solutions of the form 



J7(a:,t) = + coth[n't + f3]j X 



V 



-,(coth[l^'t + /3]-^)) ,(53) 



^wosinh[r2't + V^o sinh[r2't + /?] ' 

where ¥{'!], t) satisfies the standard BE ([5]). 

a. Shock and multi-shock solitary type solutions for the forced BE p6p can be obtained from 
the general expression 



r3(t) rsit) 

where 



aiexp[pi{x,t)] +02 ex.p[p2{x,t)] + ... + akexp[pk{x,t)] 



exp[pi(x,t)] + exp[p2(a;,i)] + ■•• + cxp[pk{x,t)] 



(54) 



Pi{x,t) = -a,;^^^a; + yT-3(0 +P°, i = l,2,...,fc, 
and r3{t) is given by ((^ . r3(t) is given by ((5^ . 

b. Generalized triangular wave solution of BE ([55]) is of the form with r(t) = ?'3(i) and 
r(t) =T3(t), t > 0. 

c. Generalized N-wave solution of BE ([36]) for the over damping case can be written using 
formula p7| . 

d. The forced BE with initial conditions [/m(x,0) = —{m/x), m = 1,2,3..., has rational 
type solutions of the form 

= If .-e-,lfi) ( !!^^M^) , „ . ,,3,.... 

where r3(t) and r3(t) are given by ([5T|) and ([5^ . respectively. 

From the general form of the solution (1551) . clearly if the function y(77(x, t), T(t)), is bounded 
for i ^ 1, then the long-time behavior of U{x,t) is described by limt„j.oo U{x,t) — (—7/2 + il')x, 
and the limiting function is a stationary solution of the BE p6p . 
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3.2.4 Variable coefficient case with cjq < 0. 

In the study of damped harmonic oscillator, usually 7 and ujq are positive parameters leading to 
the critical, under and the over damping cases, which we already discussed. However, we can take 
also ojq — — w^, a) > 0, 7 > 0, and consider the Burgers equation 

Ut+jU + UUx = ^e^^*Uxx+^^'^x, -00 < a; < 00, t > 0. (55) 

In that case, we have the IVP r + jr — uj^r = 0, r(0) = tq 7^ 0, r(0) = 0, which has solution 

r{t) = '"o^e"^ cosh[f2i + /3], 

where = 1/0)2 (72/4)^ and /3 = tanh~^(7/2i7). Then, BE 1^ has solutions of the form, 

U{x,t) = (-(7/2) + fitanh(m + /3))x 

( '-:"\] V ( ? ( -^'f ) , 4 f .a„M& + ffl - V (56) 

wVcosh(m + ^)/ Vcosh(m + /3)/ \^ ' 2f]yy 

In that case, when 7 — )■ 0, one has 17 — a), and /3 — > 0. If T^(?7, r) is as in the previous cases, then 
one can see that the solution (1551) of the variable parametric BE (|55p approaches the solution (|5T|) 
of the constant coefficient forced BE pop . which in turn approaches the standard BE, when w — > 0. 

4 Conclusion 

Nonlinear PDE's are important tool to study many physical and natural phenomena. The reach 
and complicated structure of their solutions corresponds to the reach and complicated character 
of the real world systems. In general, it is impossible to solve analytically a given nonlinear PDE. 
Only special class of integrable models admit exact solutions and these models play crucial role 
on revealing the nature of many physical phenomena, as well as provide convenient schemes to 
develop perturbation theory and test numerical methods. 

In this article, we introduced exactly solvable nonhomogeneous Burgers equation with specific 
time variable coefficients, and analytic solution of the initial value problem was obtained in terms 
of a corresponding second order linear differential equation. Burgers equations with constant 
damping, exponentially decaying diffusion coefficient and a forcing term linear in space variable 
were studied as particular cases. Generalized shock solitary waves, triangular waves, A'^-waves 
and rational type solutions were found explicitly and graphically illustrated. Special properties 
like fusion of shocks in traveling wave solutions, and motion of poles of rational type solutions 
were observed. In addition, we shortly discussed the limiting case of the parametric equations, 
and the long-time behavior of their solutions. We remark that, our results can be applied also to 
study a wide class of variable parametric Burgers and KPZ equations related with the classical 
Sturm-Liouville problems for the orthogonal polynomials, |27j . 
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Finally, we note that, there are different approaches to study the variable parametric Burgers 
problems posed in this article. The one, which we used here, is transforming the nonhomogeneous 
Burgers equation with variable coefficients to a standard Burgers equation, and then applying Cole- 
Hopf linearization. Another approach is a direct linearization of the variable parametric Burgers 
equation in the form of a variable parametric heat equation, which in turn can be transformed to 
a standard heat equation or can be solved using the evolution operator method. These problems 
are under consideration. 
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